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We present results of studies of charge disproportionation (CD) and spin differentiation in 
Na^OoC^ using the correlated band theory approach (local-density approximation+Hubbard 
U:LDA+U). The simultaneous CD and gap opening can be followed through a first order charge 
disproportionation transition. By comparison with experiments, we propose a value of the Coulomb 
repulsion strength XJ(x) that has significant dependence on the carrier concentration x, for which we 
obtain a phase diagram. The connection between the oxygen height and effects of on-site correlation 
is also reported. 

PACS numbers: 71.20.-b,71.20.Be,71.27.+a 



Since the discovery of superconductivity with T c 
~ 4 K in water-intercalated Na^OoC^ (NxCO) by 
Takada et al. in the early of 2003,0 several 
hundred papers have been published. The sys- 
tem shows several similarities to high T c supercon- 
ducting cuprates: layered transition metal oxide; 
vicinity of (presumed) Mott insulator C0O2; varia- 
tion with carrier concentration x; a "superconduct- 
ing dome." However, it also has differences from 
cuprates: electron-doped system from Mott insula- 
tor phase, actually nonmagnetic metallic sister x = 
phase, triangular lattice which can be frustrated, dif- 
ferent superconducting dome shape, and low T c . 

In addition to its interesting properties in view of 
superconductivity, its normal state has a rich phase 
diagram depending on carrier concentration x. 2] For 
x < 0.5 the Pauli-like susceptibility has been ob- 
served, while for x > 0.5 the susceptibility shows lo- 
cal moment (Curie- Weiss) behavior. For 0.75 < x < 
0.85, NxCO shows antiferromagnetic (AFM) align- 
ment of ferromagnetically ordered CoOg layers. 
Specifically at x=0.5, the system has been observed 
to undergo a charge disproportionation (2Co 3 5+ — * 
Co 3+ +Co 4+ ) and metal-insulator transition at 50 K, 
while the rest of the phase diagram is metallic. 

We have focused on normal state electronic and 
magnetic structure of NxCO, with specific interest 
in ordering near the ground state, and character- 
izing the strength and consequences of correlation 
effects. In transition metal oxides in general, the 
mechanism of charge ordering or charge dispropor- 
tionation (CD) is an important issue, and our work 
provides new inroads in the understanding of the 
mechanisms of disproportionation. 

Our calculations were based on the supercell 
approach, |j, |5| using the full-potential local-orbital 
method (FPLO).[g Both popular schemes for 
LDA+U functional were monitored. (The intra- 
atomic exchange integral J=l eV was left un- 
changed.) Both have the same Hubbard-like density- 



density interaction, but differ in just how to subtract 
out "double counting" of the pair interaction. One 
way, the so-called Fully Localized Limit (or Atomic 
Limit) ,0 is an atomic-like treatment and is appro- 
priate for large on-site Coulomb repulsion U, e.g. 
strong localization of the correlated orbital. Another 
is often called the Around Mean Field Q and is more 
appropriate when U is not so strong. In this sodium 
cobaltate, our calculations have shown similar re- 
sults from both schemes. 

Within the local density approximation (LDA), 
all Co magnetic moments are nearly equivalent, that 
is, there are only Co t - 4 ~ x ) + ions. The distorted Co- 
06 edge-sharing and the layered structures make 
t% g manifold split into & g +e' g . Upon introducing 
U (there is no agreement on the best values), both 
magnetic moments decrease linearly at small U. At 
a critical value U c , the moments change discon- 
tinuously due to a first-order transition: the state 
with similar magnetic moments no long is a so- 
lution, and the new state shows changes in mo- 
ments just as expected from charge disproportion- 
ation Co( 4 - :c )+^.tCo 3 + + (1 - x)Co 4+ . This change 
occurs simultaneously with metal-insulator transi- 
tion associated with the Mott transition of the a g 
orbital. Before CD, the a g band has a similar band 
center and width as the doubly degenerate e' g band. 

One example of this CD is shown in Fig. for 
x = 0.5, note that it is possible to follow hystere- 
sis within the region 3.2 < U < 3.6 eV and outside 
this range only one of the states, CD or undispropor- 
tionated (UD), is obtained. The two solutions can 
be observed in another way, using fixed spin moment 
(FSM) calculations within LDA+U in the critical re- 
gion (bottom panel of Fig. In this method, the 
energy is calculated versus a constrained value of the 
total moment, and the energy versus total magnetic 
moment plot becomes two curves, in which one state 
or the other is reached depending on starting point. 
Compared with LDA FSM results by Singh, [|| our 
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FIG. 1: Upper panel: Effect of the on-site Coulomb re- 
pulsion U on the magnetic moments for ferromagnetic 
ordering in Nao.sCoO^. The change shows a first-order 
transition (and accompanying hysteresis) in the critical 
region (3.2 eV to 3.6 eV). Metal- insulator transition si- 
multaneously occurs with CD (Col^Co 3+ , Co2-^Co 4+ ) 
due to nature of the first-order transition. Dashed and 
solid lines indicate the CD and UD states, respectively. 
Bottom panel: Fixed spin moment calculation of energy 
vs. total magnetic moment M (in ^is/2Co) at U=3.3 eV 
in Nao.5Co02- The left and right sides of the y-axis de- 
note the total energy difference and difference between 
minority and majority Fermi energy (effective applied 
field), respectively. The arrows mark zeroes of the Fermi 
energy difference, which correspond to extrema of the 
energy. 



LDA+U FSM shows two differences that are worthy 
of note. First, at M — 0/j,b, this system is antifer- 
romagnetic (AFM) with Co magnetic moment 0.34 
fi b ; Singh started from a nonmagnetic state at M=0 
and kept all Co ions identical. Secondly, at small M, 
the energy vs. total moment plot is nearly flat, but 
having two extrema [stable at M = 0.16/is and un- 
stable at M = 0.5/ib}- Due to the miniscule energy 
barrier, this small net moment minimum will not be 
accessible experimentally. 
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FIG. 2: Effect of Na concentration x on the critical on- 
site Coulomb repulsion U c for CD. Allowing antiferro- 
magnetic ordering, U c drops largely by about 2 eV for 
x > 0, while the change is negligible for x — 0. The 
arrows denote drop of the effect on-site Coulomb repul- 
sion U e ff due to three-band nature for x < 0.5 (see text). 
Considering ferromagnetic ordering, comparison with ex- 
periments suggests a change of U with sharp jump at 
x — 0.5, depicted by the shaded region. The data arc 
from Ref. II II El. 



As shown in Fig. [5] and described in the caption, 
the critical value U c for disproportionation depends 
on the carrier concentration x and the type of mag- 
netic order. In the lower half of the range x < 0.5, 
U c increases linearly for ferromagnetic (FM) order- 
ing. However, allowing AFM ordering, U c (ar) shows 
somewhat complicated behavior. While U c for both 
FM and AFM is almost identical at x = 0, the value 
for presumed AFM order decreases sharply for x^O. 
From Fig. U c for AFM is smaller by 2 eV than 
that of FM for x=l/3 and £=1/2. (note U c ~ 1 eV 
for AFM at x =1/3). [The supercells we used do not 
allow for combined disproportionation and AFM or- 
dering for x =2/3.] 

Furthermore, comparison with experiments sug- 
gests that U(x) is significantly dependent on the car- 
rier concentration x. First, no observation of CD nor 
any significant correlation effects at all for x < 0.5 
indicates that U(x < 0.5) is less than U c . How- 
ever, the study of Gunnarsson and coworkers shows 
that in a N-multiband system the effective Coulomb 
repulsion becomes U e ^ = U / \/N. [lll| Below U c , the 
a g band is only distinguishable from the t2 g mani- 
fold (in energy) by lying slightly higher, and hosting 
most of the 1 — x holes. As a result the system at 
x < 0.5 maintains multiband (three-band) behav- 
ior, resulting in U^ // =U c /\/3- 1.4 eV and 1.7 eV 
for x = and x — 1/3, respectively. Since only a 
tiny band gap (^15 meV) is observed at x = 0.5.|12j 
\J(x = 0.5) must be near the minimum value for CD, 



2 



* 




- * 




* 
* 


I 

\ 

v 


* 

- # 


* x + — 

# \ * + * 


. Metal X Half-Metal 


* 


* * + + 

"N. 







0.2 0.21 0.22 0.23 0.24 
Oxyen height (z ) 

FIG. 3: Phase diagram depending on oxygen height (zo) 
and U for ferromagnetic ordering at x — 0. (Energeti- 
cally, zo = 0.235 is favored.) Lattice constant c = 4.2509 
A is used. The symbols describe metallic (*), half- 
metallic (+), and insulating (diamond) states. The data 
are from Ref. ||. 

i.e. 3.5 — 4 eV. For x > 0.5, the observed Curie- Weiss 
susceptibility indicates (1 — x)Co 4+ and ccCo 3+ ions, 
so that the charge is disprop orti onated even though 
the system is a good metal. O Thus, U(x > 0.5) 
is above 3 eV, perhaps over 4 eV. We will discuss 
this £=0.5 state elsewhere. It remains for further 
research to decide how accurately LDA+U should 
be for the band gap, or whether dynamics of the 
correlations are important. 

Another issue which is drawing attention in these 
cobaltates is whether the change of O height (i.e., 
Co-0 bond-length) is important as the carrier con- 
centration x increases. Johannes et al. have calcu- 
lated the oxygen height zq by energy minimization, 
and concluded that the ch ang e has little effect in 
the range 0.3 < x < 0.75. |14j However, their cal- 
culations, and therefore conclusions, are confined to 



the local density approximation results and predic- 
tions. We have carried out a study at x=0 to address 
whether correlation effects described by the LDA+U 
method are affected by the oxygen height. Figure |3 
shows phase diagram in the relevant portion of the 
U - zo plane in the relevant range of U and allowing 
FM ordering of the system. 

For z > 0.21, the system undergoes a metal- 
insulator transition, implying CD, as U is increased. 
The critical value U c depends little on O height. In 
fact, AU C is only 0.1 eV between z$ = 0.22 and 
zq = 0.23, which is the same variation of O height 
as measured between x = 0.3 and x = 0.7. This re- 
sult suggests that the critical value U c for CD is lit- 
tle affected by O height in this system. In addition, 
it is worthwhile to note that a metal to half-metal 
transition occurs even for very small U before the 
transition to insulator (for example, at U=l eV for 
z = 0.23). 

Note however that for a height xo=0.20, no CD 
is obtained for U up to 5.5 eV, much different from 
the behavior for zo > 0.21. So while this bondlcngth 
lies outside the accessible range, it does show that 
there is a regime where correlation effects depend 
very strongly on the Co-0 bond length, and that the 
system is not so far from that regime. Finally, we 
note that the transitions with varying U discussed 
here should be of special interest for high pressure re- 
search because, for a system near U c , varying U will 
be analogous to applying pressure to change U/W 
ratio (W is the correlated orbital's bandwidth). 
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